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Abstract 

The eigenvalues of a non-Hermitian Hamilton operator are complex and provide not only the 
energies but also the lifetimes of the states of the system. They show a non-analytical behavior 
at singular (exceptional) points (EPs). The eigenfunctions are biorthogonal, in contrast to the 
orthogonal eigenfunctions of a Hermitian operator. A quantitative measure for the ratio between 
biorthogonality and orthogonality is the phase rigidity of the wavefunctions. At and near an EP, 
the phase rigidity takes its minimum value. The lifetimes of two nearby eigenstates of a quantum 
system bifurcate under the influence of an EP. When the parameters are tuned to the point of 
maximum width bifurcation, the phase rigidity suddenly increases up to its maximum value. This 
means that the eigenfunctions become almost orthogonal at this point. This unexpected result is 
very robust as shown by numerical results for different classes of systems. Physically, it causes an 
irreversible stabilization of the system by creating local structures that can be described well by a 
Hermitian Hamilton operator. Interesting non-trivial features of open quantum systems appear in 
the parameter range in which a clustering of EPs causes a dynamical phase transition. 
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I. INTRODUCTION 


In its simplest form, the wavefunction of a quantum many-particle fermionic system 
can be written as a Slater determinant that satishes the anti-symmetry requirements and 
consequently the Pauli principle. However, some corrections in the wavefunctions are usually 
necessary which are caused by their possible mixing. In nuclei, for example, corrections arise 
from the residual interaction between the particles of the system, with the consequence that 
the wavefunctions of all the states are strongly mixed. Another source is the interaction 
between system and environment into which the system is embedded. In this case, the 
states of the system interact via the common environment due to which their wavefunctions 
may be modihed. 

The natural environment of a localized quantum mechanical system is the extended con¬ 
tinuum of scattering wavefunctions in which the system is embedded. This environment can 
be changed by means of external forces, however it can never be deleted. It causes some 


communication between distant levels, for details see the recent review [l| in which theoret¬ 
ical results are confronted with experimental results. The theoretical results are obtained 
by using a non-Hermitian Hamilton operator "H for the description of the system, which 


contains explicitly the interaction between system and environment [^. The eigenvalues 
and eigenfunctions of 'H differ essentially from those of a Hermitian Hamilton operator: the 
eigenfunctions are biorthogonal and the eigenvalues may show, as function of a certain pa¬ 
rameter, deviations from Fermi’s golden rule. The differences between the eigenvalues and 
eigenfunctions of H and those of a Hermitian operator H appear, above all, at and near to 
singular points at which two eigenvalues of T-L coalesce and the corresponding eigenfunctions 
differ from one another only by a phase. These singular points do not have any equivalent 
in the mathematics based on Hermitian Hamiltonians since the different eigenfunctions of 
H are orthogonal (in contrast to those of which are biorthogonal). These singular points 
are called usually exceptional points (EPs), according to j^. 

The role of EPs in physical systems is considered in many papers during last years, e.g. 
4-8|, see also the review and the book 9|. In the present paper, we are interested in the 
relation of these singular points to phase transitions occurring in open quantum systems. 
Such a relation is hrst discussed theoretically some years ago in 


lol l and 


11| . however 


without rigorous consideration of the biorthogonality of the eigenfunctions of "H. The same 
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holds true for the papers on the dynamical phase transition observed experimentally 
and theoretically in the spin swapping operation in atomic systems. Only recently, the 
mixing of the wavefunctions at an EP and its relation to a dynamical phase transition is 


studied for an open quantum system with more than two states 


12l |. The states at both 


sides of the phase transition are not analytically related to one another, meaning that this 
characteristic feature of any phase transition is fulhlled. The theoretical results are compared 


fl. 


to experimental results in the recent review 

In contrast to these differences between Hermitian and non-Hermitian quantum physics 
at and near to an EP, the Hamilton operator of the Schrodinger equation of an open quantum 
system is almost Hermitian, to a good approximation, far from an EP. Here the eigenfunc¬ 


tions of the non-Hermitian Hamilton operator are almost orthogonal 


13| : and the system 


can be described quite well by a standard Hermitian Hamilton operator. We underline how¬ 
ever that the Hamiltonian remains dehnitely non-Hermitian since the non-Hermiticity arises 
solely from the fact that the function space of the localized part of the system is a subsystem 
of the total function space: the localized part, we are interested in, is embedded into an 
extended environment of scattering wavefunctions, see Fig. 1 in Ij. The eigenfunctions of 

fl fl 

a non-Hermitian Hamilton operator are always biorthogonal [2| (see also Sect. 3.2 in [![). 
This includes the case with well-separated resonance states as shown theoretically [l^ and 


experimentally 


id- 


It is the aim of the present paper to consider the influence of EPs onto the eigenfunc¬ 
tions (including their phases) of a non-Hermitian Hamiltonian in detail. We are interested, 
above all, in the behavior of the eigenfunctions and their phases in approaching an EP and, 
furthermore, in the role they play in a dynamical phase transition. 

First, we sketch in Sect. m the characteristic features of the eigenvalues and eigenfunc¬ 
tions of a non-Hermitian 2x2 Hamiltonian which are known in literature, and provide the 
dehnitions of typical values such as, among others, the phase rigidity. This is a quantitative 
expression for the biorthogonality of the wavefunction and thus, ultimately, of the degree of 
opening of the system. 

In the following Sect. IHlI we consider the conditions for the appearance of EPs in different 
systems consisting of two states. Analytical results are obtained and discussed for cases with 
two EPs. Then, we consider in Sect. llVI the more general case of a system with more than 
two states where a clustering of EPs may occur. In Sect. |V]we provide numerical results for 
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systems under typical conditions and compare them with analytical results, if such results 
exist. We discuss the obtained results in Sect. m Some of the results are expected and 
agree with our general understanding of open quantum systems. Other results are completely 
unexpected. These results are considered and discussed in detail. They allow us to receive 
a deeper understanding of dynamical phase transitions. We conclude the paper with some 
remarks on the stabilization of open quantum systems due to the existence of EPs; and 
discuss the possibility to describe them by means of a Hermitian Hamiltonian. 


II. EIGENVALUES AND EIGENFUNCTIONS OF THE NON-HERMITIAN 
HAMILTONIAN ^(2) 

In order to study the interaction of two states via the common environment it is most 
convenient to start from the symmetric 2x2 non-Hermitian matrix 


m 


^( 2 ) _ 

with 7 i < 0 for decaying states 


^ = Cl + |7i 


V 


1^21 



( 1 ) 


16| . The Ui 2 = 0021 = oo stand for the coupling matrix 
elements of the two states via the common environment which are, generally, complex j^. 
The diagonal elements e* of ([T]) contain the energies e* and decay widths 7 * of the two 
states when Uij = 0 , i.e. they are the two complex eigenvalues Si (z = 1 , 2 ) of the non- 
Hermitian operator "Hq that describes the system without any coupling of its states via the 
environment. In the present paper, we take 0212 = 0021 = 00 , and the self-energy of the states 
is assumed to be included into ei and £ 2 - We underline here that the Hamiltonian is 
completely non-Hermitian (see also [^) in difference to the many non-Hermitian operators 
used in the literature for the description of open quantum systems. These operators consist 
mostly of a Hermitian part to which a non-Hermitian part is added as a perturbation. 

The model ([T]) seems to be very simple. This is however not true from a mathematical 
point of view. The point is that singularities are involved in the model (the so-called EPs) 
which are known in mathematics for many year^ see [3]. They are considered in physics only 
recently. They cause counterintuitive results |2| which seem to be, at hrst glance, wrong. 
We will discuss them in the following. 

In the case a; = 0, the real part Ek of the eigenvalues of does not differ from the 
original energies e^. It follows, under this condition, Ei ^2 = ^ (61 + 62 ) ± | (61 — 62 ) = 61 ^ 2 - 
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A corresponding relation holds for the relative to the original 7 ^. 

Most interesting properties of are the crossing points of two eigenvalue trajectories. 
Since here the two states coalesce at one point, the influence of all the other states of the 
system on the interaction of these two states can be neglected. Therefore, ([T]) describes the 
characteristics of open quantum systems that may be related to these points, in spite of its 
small rank. 

The eigenvalues of are, generally, complex and may be expressed as 


^1,2 = ± Z ; 


^ = 2 V “^ 2 )^ + 4a;2 


( 2 ) 


where Ei and T* stand for the energy and width, respectively, of the eigenstate i. Also here 
Tj < 0 for decaying states 1^. The two states may repel each other in accordance with 
Re(Z), or they may undergo width bifurcation in accordance with Im(Z). When Z = 0 the 
two states cross each other at a point that is called usually exceptional point (EP) j^. The 
EP is a singular point (branch point) in the complex plane where the S-matrix has a double 
pole j^. 

The eigenfunctions of any non-Hermitian operator "H must fulhll the conditions = 

£i\^i) and (TjI'H = £i{^i\ where £i is an eigenvalue of 1-i and the vectors |<hj) and (Tjl de¬ 


note its right and left eigenfunctions, respectively j^. When is a Hermitian operator, the 
£i are real, and we arrive at the well-known relation (Tjl = (<hj|. In this case, the eigenfunc¬ 
tions can be normalized by using the expression For the symmetric non-Hermitian 

Hamiltonian however, we have (Tjl = ($*]. This means, that the eigenfunctions are 
biorthogonal and have to be normalized by means of (<h*|<hj). This is, generally, a complex 
value, in contrast to the real value (<hi|<hj) of the Hermitian case. To smoothly describe 
the transition from a closed system with discrete states, to a weakly open one with narrow 
resonance states, we normalize the according to 


{•S>-|4>,> = % 

(for details see Sects. 2.2 and 2.3 of j^). It follows 

= Re (($.1$.)) ; A = ($*|4>*) > 1 

and 

= i Im ((<hi|<hjyi)) = 

> 0 . 


( 3 ) 


( 4 ) 


( 5 ) 
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At the EPs, not only the eigenvalues of two states coalesce but also the two corresponding 
eigenfunctions of the non-Hermitian Hamilton operator are the same, up to a phase, 




( 6 ) 


These relations follow from analytica 


see Appendix of 


17|, Sect. 2.5 of 


as well as from numerical and experimental studies. 


and Figs. 4 and 5 in 


18| . We underline here that 


the coalescence of the two eigenvalues of a non-Hermitian operator at an EP should not 
be confused with the well-known fact that two eigenstates of a Hermitian operator may 
be degenerate. The difference consists in the relation (E]) between the two corresponding 
eigenfunctions which does, of course, not hold for degenerate states. Furthermore, the 
eigenfunctions of a non-Hermitian operator are biorthogonal while those of degenerate states 
are orthogonal. 


According to ([6]), the wavefunction $i of the state 1 jumps, at the EP, to $2 


This mathematical behavior of the eigenfunctions <hj at the singular EPs causes the main 
differences between the physics of Hermitian and non-Hermitian quantum systems. At an 
EP, 00 , \Bl\ 00 , and the influence of the environment onto the system is extremely 

large |2|. 

In (|3]), the complex value (<I>*|<I>j) is normalized to the real value 6ij with the consequence 
that the relative phase between the biorthogonal eigenfunctions of two neighbored states 
changes in such a manner that always Im(<h*|<I)j) = 0. A quantitative measure of this 
change is the so-called phase rigidity 


^ (4>*|«I>fc) 
" {^k\^k) 


= 


(7) 


of the state k which is dehned by the ratio between biorthogonality and orthogonality of 
the wavefunctions <hfc. For Hermitian systems for which (4)^1$^) = (<hfc|<hfc) holds, the 
phase rigidity is equal to unity. This is an expression of the fact that the eigenfunctions 
of Hermitian operators are orthogonal. For weakly decaying systems, where one has well- 
separated resonance states, the wavefunctions are almost orthogonal, i.e. the degree of 
biorthogonality is small. Under such conditions, Hermitian quantum physics represents a 
reasonable approximation to the description of the open quantum system. However, the 
wavefunction $1 jumps at the EP to ±i $2 and vice versa; and the phase rigidity does not 
vary continuously at the EP also in this case. 
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The Schrodinger equation with the non-Hermitian Hamilton operator is equivalent 

( 2 ) 

to a Schrodinger equation with "Hq source term 

/ 


!$,) = 


0 oJi 


14,) = »y|4,). 


( 8 ) 


y ^ji d 

This equation relates to in a non-trivial manner due to the source term. That means, 
two states i and j ^ i are coupled via the common environment of scattering wavefunctions 
into which the system is embedded. It is = Uji = ca, and the coupling between the states 
i and j ^ i vanishes when ca —?■ 0. The Schrodinger equation ([H]) with source term can be 
rewritten in the following manner 20|, 


{nf -Si) |<i>,) = ^ . (9) 

fc=l,2 m,=l,2 

According to the biorthogonality relations (jlj) and ([5]) of the eigenfunctions of ([9]) is 
a nonlinear equation, since (<hfc| <!>„,) 7 ^ 1 for k = m and (<I>fc|<hm) 7^ 0 for k ^ m. Most 
important part of the nonlinear contributions is contained in 


(Wf - En) |4„) = (4,.|H^|4„) |4„|" |4„) . 


( 10 ) 


The nonlinear source term vanishes far from an EP where (<hfc|<hfc) approaches 1 and 

approaches zero. This follows from the normalization ([3]) which 
differs only a little from the standard normalization = 1 and = 0 for a 

Hermitian Hamilton operator. Thus, the Schrodinger equation with source term is (almost) 
linear far from an EP, as usually assumed. It is however nonlinear in the neighborhood of 
an EP. 

The nonlinear terms in (iTOll cause, among others, a mixing of the wavefunctions which 
can be expressed by 




N 


E'>««?■ 


( 11 ) 


where the eigenfunctions <I>fc of are represented in the set of basic wavefunctions of 
the operator T-Ln the non-diagonal matrix elements of which vanish. For some numerical 

n 

results see [13|. This mixing of the wavefunctions appears additionally to other possible 
sources of mixing caused for some other reasons. 
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The eigenfunctions <hj and the eigenvalues Si of contain global features that are 
caused by many-body forces induced by the coupling Uik of the states i and k ^ i via the 
environment. The environment is the continuum of scattering wavefunctions and has an 
inhnite number of degrees of freedom. 


III. EXCEPTIONAL POINTS 


We consider now the behavior that arises when the parametrical detuning of the two 
eigenstates of is varied, bringing them towards coalescence. According to ([ 2 ]), the 
condition for coalescence reads 

Z = ^^{ei- 62^ - ^(71 - 72)^ + i{ei - 62)(71 - 72) + 4 a ;2 = 0 . ( 12 ) 

It follows that two interacting discrete states (with 71 = 72 = 0 and ei 7 ^ 62) avoid always 
crossing since u = uq and £1 — £2 are real in this case and the condition Z = 0 can not be 
fulhlled, 

(ci — 62)^ T > 0 . ( 13 ) 


In this case, the EP can be found only by analytical continuation into the continuum. 
This situation is called usually avoided crossing of discrete states. It holds also for narrow 
resonance states if Z = 0 cannot be fulhlled due to the small widths of the two states. The 
physical meaning of this result is very well known since many years: the avoided crossing of 
two discrete states at a certain critical parameter value 2 l|] means that the two states are 
exchanged at this point, including their populations {population transfer). 

When 7 i = 72 , and a; = icuo is imaginary^ it follows from (IT^ 


(ei — 62 )^ — 4 ojq — 0 —y ei — 62 — i 2 uq (14) 

such that two EPs appear. It furthermore holds that 

(ci — 62)^ ^ dcug — y Z ^ (15) 

(ci — 62 )^ ^ dcug —y Z E ^ (16) 

independent of the parameter dependence of ei, 2 - 
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Also in the case when the widths 71,2 are parameter dependent, ei = 62 and oj is real, we 
have two EPs. Instead of flT^ to flT6|l we have 


and 


(71 - 72)^ - 16 


0 7i — 72 = ±4a; 


(17) 


(71 — 72)^ > ISo;^ —t 2’ G O' 
(71 - 72)^ < 16 ci;^ -)■ Z G 3 ? . 


(18) 

(19) 


Eqs. da and (la. respectively, describe the behavior away from the EPs, where the 
eigenvalues £k only differ from the original ones through a contribution to the energy and 
width, respectively. The widths (or energies), in contrast, remain unchanged, and this 
situation therefore corresponds to that of level repulsion (or width bifurcation). Eqs. flTB]l 
and flT^ . in contrast, are relevant over the range between the two EPs and indicate that the 
resonance states undergo width bifurcation (or level repulsion) according to lm{Z) 7^ 0 (and 
Re(Z) 7^ 0, respectively). The bifurcation (or level repulsion) starts in the neighborhood 
of either one of the EPs, and grows to reach a maximum value at the midpoint between 
them (even though uq and u, respectively, remain hxed). The condition for maximum 
width bifurcation (or level repulsion) is fulhlled at the crossing point ei = 62 (and 71 =72, 
respectively). Physically, the bifurcation implies that different time scales may appear in 
the system, while the states are nearby to one another in energy (for details see 17|). In 


an analogous manner, level repulsion of states with similar lifetimes causes a separation of 
the states in energy. For an illustration of these analytical results see the numerical results 
presented in Sect. IVA1 Fig. [Hand Fig. [3] left panel. 

Under more realistic conditions, uj is complex, and simple analytical results like (ITTll to 
flT9|) cannot be obtained. For this case, we will provide some results of numerical studies 
in Sect. IV Al In order to understand the meaning of these numerical results, the analytical 
relations (na to f|19p and their representation in Figs. [1] and Fig. [3] left panel are very 
helpful. 

In any case, the parametric dependence of the eigenvalues Sk is non-analytical in the 
vicinity of an EP, with the widths Pfc, in particular, showing variations that are inconsistent 
with the predictions of Fermi’s golden rule (according to which the widths should increase 
with increasing coupling strength of the system to the environment; for details see j^). 
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At these points, the influence of the environment onto the system properties is extremely 
strong. In our case, the environment is the continuum of scattering wavefunctions which 
gives to the eigenstates of a finite lifetime. 

It follows from the normalization condition (|3]) that —)■ oo such that —)■ 0 

when an EP is approached. In other words, the relative phase of the two eigenfunctions 
changes dramatically when the crossing point is approached. Most significantly, as under¬ 
stood from analytical studies, as well as from numerics and experiment (see |2|, [iTI, ll9|) is, 
that, in the vicinity of the EP, the eigenfunctions differ from one another by only a phase, 
see (E]). In a recent theoretical study on a microwave cavity [^, the relations (E]) could 
be confirmed by the observation that the real and imaginary components of two nearby 
eigenstates are “swapped”, under the influence of an EP, in complete agreement with (E]). 
The non-rigidity of the phases follows, of course, directly from the fact that (<h^|<hfc) is a 
complex number (in difference to the norm ($^1$^), which is a real number) so that the nor¬ 
malization condition (E]) can be fulfilled only by the additional requirement Im(<h^|<hfc) = 0 
(corresponding to a rotation away from the complex plane). Here, the two different states 
of the system develop, according to (El), a coupling through the continuum, a quantitative 
measure of which is the phase rigidity. Thus, the biorthogonality of the eigenfunctions 
causes perceptible physical effects in the neighborhood of an EP. 

Generally speaking, the phase rigidity takes values between zero and one, with the value 
Tfc = 1 for Hermitian systems. Near to an EP in a non-Hermitian system, however, the 
two eigenfunctions differ from one another only by a phase, according to (E]), so that r*, -C 
1. This non-rigidity of the eigenfunction phases is the most important difference between 
Hermitian and non-Hermitian eigenfunctions. Its meaning cannot be overestimated. On the 
one hand, the lack of phase rigidity near to an EP leads very naturally to the appearance 
of nonlinear effects in the Schrodinger equation (fTOj) with source term which describes an 
open quantum system. On the other hand, the impact of the environment on the (localized) 
system is extremely strong at the EP. Since the environment is the continuum of scattering 
wavefunctions with an infinite number of degrees of freedom, this impact may induce phase 


transitions as discussed in 


12 |. 
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IV. CLUSTERING OF EXCEPTIONAL POINTS 


According to mathematical studies, more than two eigenvalues of a non-Hermitian oper¬ 
ator 'H may cross in one point, the so-called higher-order EP. This crossing point is however 
a point in the continuum and therefore of measure zero |^. In this respect, it does not 
differ from the second-order EP which is the crossing point of two eigenvalues considered 
in the foregoing Sect. [TTl That means, a higher-order EP can not directly be identihed in 
a realistic physical system. Nevertheless, it influences the dynamics of an open quantum 
system in a similar manner as a second-order EP does; see the discussion in the foregoing 


Sect. [IT1_ 


In 


or a two-level system. 


12l | , the influence of a third state onto the two eigenvalues and eigenfunctions of a non- 


Hermitian Hamilton operator that cross at an EP, is investigated. As a result, more than 
two states of a realistic physical system are unable to coalesce at one point since, in a certain 
hnite parameter range around the original second-order EP, the wavefunctions of the two 
states are mixed. When the third state approaches this parameter range, it crosses or avoids 
crossing therefore with states that differ from the original two states. Accordingly, new 
EPs appear and the areas of influence of different EPs overlap. Altogether, the different 
EPs amplify, collectively, their impact onto physical values; and the wavefunctions of all 
states are strongly mixed in the basic wavefunctions of "Hq- This effect is nothing but 
some clustering of EPs, wherewith the characteristic fact is expressed that the ranges of 
the influence of different second-order EPs overlap in a hnite parameter range around a 
higher-order EP. 

In the following, we will study the mixing of the wavefunctions and, above all, the phase 
rigidity dehned in, respectively, CD and ([7]), in the case of clustering of EPs. To this aim 
we consider the non-Hermitian Hamiltonian 
/ 


Jlin) = 


ei = Cl -f i^ 
(^21 


V 


^nl 


^12 

€2 = e2+ i Y 
0 
0 


0 


0 


^1% 

0 

0 


f = P A- 1 




( 20 ) 


/ 


with n = 3 or 4 nearby states coupled to one common continuum (the hrst channel). As 
in (II]), the e* and denote the energies and widths, respectively, of the n states without 


account of the interaction of the different states via the environment. The Uij = Uji simulate 
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the interaction of the two states i and j via the common environment. In the simulation 
fl20|) . we used the doorway concept used in nuclear physics: the n states with the decay 
widths 7j/2 can be simulated by one doorway state with large decay width 71 and n — 1 
states with small (almost vanishing) decay widths 7*^1. Then (according to the doorway 
concept), the doorway state is coupled to both the environment and the remaining n — 1 
states, while the remaining states are coupled to the environment only via the doorway 
state (due to their small decay widths and the fact that they are distant from EPs). The 
coupling strength uj between system and environment is not varied in our calculations, and 


22 |. 


the number of parameters for the widths and energies of all n states is 2n 

The eigenvalues of fl2(l can be obtained in analogy to ([2]). The eigenfunctions are 
biorthogonah We normalize them according to ([3]). Further (E]) holds at an EP. The values 
Ai and \Bl\ dehned in (jl]) and (jS]), respectively, express how near the system is to an EP 
at the considered parameter value. In the numerical calculations, these values can be seen 
directly by studying the mixing coefficients \bij\ defined in (1TT|1 . Also the corresponding 
phase rigidities r* of the different states can be determined numerically by using (JTj). 

Near to the different EPs, the Schrodinger equation contains nonlinear contributions 
according to flTOl) due to the source term that describes the coupling between system and 
environment. Accordingly, the whole parameter range in which a clustering of EPs occurs, 
is controlled by nonlinear contributions to the Schrodinger equation, the values of which 
vary because of their dependence on the concrete parameter value. They vanish only far 
from this regime with a clustering of EPs. 


V. NUMERICAL RESULTS 

A. N = 2 states 

In Figs. [1] to m we show the results of numerical calculations performed with the pa¬ 
rameters given in Table [B Most impressive is that all results for the phase rigidity show the 
same behavior in spite of the different parameters and the fundamental differences in the 
eigenvalue pictures. In all cases, the phase rigidity approaches the value —)■ 0 at the 
position of the EP while it approaches sharply the value —)■ 1 when width bifurcation and 
level repulsion, respectively, is maximum. These changes occur without any changes of the 
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TABLE I: The parameters used in the calculations 


Figure 

ei 

62 

7i/2 

72/2 

UJ 

Fig. [Ha-d 

2/3 

2/3 Ad 

-0.5 

-0.5 

0.05 i 

Fig. [He-h 

1/2 

1/2 

-0.5 

-0.5 a 

0.05 

Fig. [2ja-d 

2/3 

2/3+ d 

-0.5 

-0.55 

0.025 (1+i) 

Fig. [2je-h 

0.55 

0.5 

-0.5 

-0.5 a 

0.025 (1+i) 

Fig. [3la-d 

0.5 

0.5 

0.05 a 

-0.05 a 

0.05 

Fig. [3le-h 

0.55 

0.5 

0.05 a 

-0.05 a 

0.025 (1+i) 

Fig. [Ha-e 

0.5 

a 

-0.05 

-0.06 

0.05 (A+i) 

Fig. [If-j 

0.5 

0.51 

-0.5 

—0.3 a 

0.05 (1 + Ti) 


coupling strength oj between system and environment as can be seen from the parameter 
values given Tab. [B 

Results for examples with two EPs according to Eqs. flTTll to flTHl) and flTTl) to flT^ . 
respectively, are shown in Fig. [H The numerical results agree with the analytical ones: 
the phase rigidity approaches zero at the two EPs; while in between the EPs, we see width 
bifurcation in the hrst case (Fig. [T]b) and level repulsion in the second case (Fig. [He). The 
result Tfc —)■ 0 at every EP corresponds to the expectation of theory j^. However, there is 
an unexpected sharp transition to —)■ 1 at the point of maximum width bifurcation or 
maximum level repulsion. This means that here the two eigenfunctions of are (almost) 
orthogonal to one another. Far from the critical region, the phase rigidity approaches the 
value 1 according to the fact that the influence of the environment onto the system can be 
neglected, to a good approximation, far from EPs. 

The results for the more realistic case with complex coupling strength uj are shown in Fig. 
O Here, only one EP appears. In a hnite parameter distance from the EP, we see maximum 
width bifurcation (Fig. |2lb) and maximum level repulsion (Fig. Ele), respectively. Again, 
rfc —?• 0 at the EP and r*, 1 at maximum width bifurcation or maximum level repulsion. 

The results in Fig. [3] show the eigenvalue and phase rigidity pictures for the case when 
not only loss (as in Figs. [Hand[2]) appears but also gain is a possible process. Fig. [3] left 
panel shows the case with balanced loss and gain, corresponding to P^ = 0 in the hnite 
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FIG. 1: Energies Ei (full lines) (a,e), widths Ti/2 (b,f), phase rigidity ri = r 2 (c,g), and l — Vk (d,h) of the 
two eigenfunctions of as a function of the distance d between the two unperturbed energies ei and 62 - 
The parameters are given, respectively, in the first and second row of Tab. ID The dashed lines in (a,e) and 
(b,f) show, respectively, the and 7^/2 trajectories. 


parameter range between the two EPs (see Fig. Ob). Formally, this case is similar to those 
discussed recently in many papers related to non-Hermitian operators with PT symmetry 


241. The PT-symmetry 


whose eigenvalues are real in a hnite parameter range, see e.g. 23|, 
breaking is caused by EPs. 

The results shown in Fig. |3] left panel have the same characteristic features as those 
shown in Fig. [1] right panel. The same holds true when the coupling strength u is complex 
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FIG. 2: The same as Fig. [T]but with the parameters given, respectively, in the third and fourth row of 
Tab. n 

(Fig. |3] right panel as compared to Fig. |2] right panel). 

In Fig. m some resnlts are shown with, respectively, almost imaginary (left panel) and 
almost real (right panel) conpling strength uj. We see again the characteristic sharp transition 
Tfc —)■ 0 in approaching the EP and —)■ 1 at another parameter valne at which we have 
maximnm width bifnrcation and maximnm level repnlsion, respectively. Additionally, we 
show in Fig. IHthe mixing of the wavefnnctions expressed by the coefficients \hki\ which are 
defined in (fTTj) . At the EP, |6fcz| —t cxo as shown in [^. As in the other fignres with complex 
conpling strength u (Figs. |2] and [3] right panel), there is only one EP. The point of maximnm 
width bifnrcation and maximnm level repnlsion, respectively, appears at a finite parameter 
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FIG. 3: The same as Fig. [T]but with the parameters given, respectively, in the fifth and sixth row of Tab. 

in 


distance from the EP. In this parameter region, the two wavefnnctions are strongly mixed. 
The mixing remains when ^ 1 is approached. Beyond ^ 1, we see the hint to a 
nearby EP (0 < r*, < 1) which limits the extension of the total critical parameter region. 
Beyond this critical parameter region, the wavefnnctions approach their original orthogonal 
character. The physical meaning of the mixing of the wavefnnctions nnder the inflnence of 


EPs is discnssed in detail in 


12 


13| . Here, we nnderline only that the wavefnnctions are 


mixed when r*, —>■ 1, i.e. when they are almost orthogonal in the critical parameter region. 
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FIG. 4: The same as Fig. [T]but with the parameters given, respectively, in the seventh and eighth row of 
Tab. m In (e) and (j), the mixing coefficients \bij\ of the wavefunctions, defined by (fTTl) . are shown. 


B. N > 2 states 


In Fig. El we show the influence of a ’’third” state onto the eigenvalue picture and the 
phase rigidity around an EP in a system that is symmetric around the crossing point of the 
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FIG. 5: Energies Ei (full lines) (a,e), widths Vij^ (b,f), phase rigidity (c,g), and 1 — n (d,h) of the 
two eigenfunctions of (left panel) and of the three eigenfunctions of (right panel), respectively, 
as a function of the parameter a. The dashed lines in (a,e) and (b,f) show, respectively, the and 7^/2 
trajectories. The parameters are uj = 0.01 i; ei = 1 — 1/2 a; 62 = a; 63 = —1/3 + 3/2 a (e-h); 
71/2 = 72/2 = -0.495; 73/2 = -0.4853 (e-h). 

Ci trajectories. New EPs can be identified in the three-level case in the eigenvalne pictnres 
and, correspondingly, in the values of the phase rigidity (r* < 1). The two EPs near to the 
crossing point of the e, trajectories are well expressed, and we see clear hints to the existence 
of two other distant EPs. Most interesting result is the sharp transition from a reduced value 
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FIG. 6: The same as Fig. [5]but lo = 0.005(1+ i) and 71/2 = —0.5; 72/2 = —0.51; ei = 0.5; 62 = a (left 
panel) and 71/2 = —0.5; 72/2 = —0.505; 73/2 = —0.51; ei = 0.5; 62 = a; 63 = 2 a — 0.5 (right panel). 


Tj < 1 of the phase rigidity (which is characteristic of an EP in the neighborhood) back to 
Tj ^ 1 at maximnm width bifnrcation. This abrupt transition occurs in both, the two-level 
system and the three-level one. At this parameter value, the eigenfunctions of the non- 
Hermitian Hamiltonian as well as those of become almost orthogonal. At and 
around this critical parameter value, the wavefunctions of the states are strongly mixed in 
both, the two-level case and the three-level case 12|. 

Fig. O shows results for a more realistic case with complex coupling strength uj between 
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FIG. 7: Energies Ei (full lines) (a,f), widths Ti/2 (b,g), phase rigidity (c,h), 1 — (d,i), and mixing 

coefficients \bij \ of the wavefunctions, defined by (fTTl) . of the four eigenfunctions of as a function of the 
parameter a . The dashed lines in (f,g) show, respectively, the and 7^/2 trajectories. The parameters are 
uj = 0.01 i] Cl = 1 — a/2; 62 = a; 63 = —1/3 + 3/2 a; 64 = 2/3; 71/2 = 72/2 = —0.4950; 73/2 = 
—0.4853; 74/2 = —0.4950 (left panel) and w = 0.005 (1 + i)] ei = 0.5; 62 = a; 63 = 2a — 0.5; 64 = 
1-a; 71/2 = -0.5; 72/2 =-0.505; 73/2 =-0.51; 74/2 =-0.505 (right panel). 

system and environment. In the three-level case, hints to the existence of different (distant) 
EPs can be seen. The sharp transition r* ^ 1 at maximnm width bifnrcation however 
appears very clearly not only in the two-level case bnt also in the three-level case. 

The resnlts shown in the Figs. |5] and |6] characterize the generic behavior of many-level 
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open quantum systems under the condition that the width bifurcation is maximum. The 
results are conhrmed by those which we received in many other calculations performed with 
different parameters (including for systems which allow loss and gain); or with a larger 
number N of states. For illustration we show numerical results obtained for = 4 states 
in Fig. [71 Again we see hints to several EPs as well as the sharp transition of rj to a value 
almost 1 at a critical parameter value. Here, the wavefunctions of all four states are strongly 
mixed. 

VI. DISCUSSION OF THE RESULTS 

A. N = 2 states 

Most surprising result of our study is the strong parameter dependence of ru in a certain 
critical region around an EP. The variation —)■ 0 in approaching an EP is expected from 
Eq. dHD- The rapid variation ^ 1 when, respectively, the maximum width bifurcation 
and level repulsion is approached, is however unexpected. The width bifurcation or level 
repulsion starts at the EP without any enhancement of the coupling strength between system 
and environment. It is driven exclusively by the nonlinear source term of the Schrodinger 
equation fllUl) which describes the open quantum system. When —)■ 1 the wavefunctions 
of the two states, the eigenfunctions of which coalesce (up to a phase factor) at the EP 
(correspondingly to ^ 0), become almost orthogonal to one another. The transition from 
rfc —)■ 0 (at the EP) to r*, —)■ 1 (where width bifurcation and level repulsion, respectively, is 
maximum) occurs suddenly as function of the varied parameter in all our calculations. The 
wavefunctions of the two states behave smoothly, i.e. they remain mixed also when —>■ 1. 

The physical meaning of this result consists in the fact that a stabilization of the localized 
part of the system occurs when the interaction u between system and environment is strong 
enough, i.e. when it is of the same order of magnitude as the widths 7j/2. Then, in the case of 
Eqs. oa to f|16p . one of the two states receives a very short lifetime due to width bifurcation, 
and becomes almost indistinguishable from the states of the environment. Although this 
process seems to be reversible according to the hgures for the eigenvalues, this is in reality 
not the case. The processes occurring in approaching —)■ 1, take place, as mentioned 
above, by means of the nonlinear source term of the Schrodinger equation near an EP. Due 
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to these processes, the long-lived state has “lost” its short-lived partner with the consequence 
that the two original states cannot be reproduced. This evolution is therefore irreversible. 
The long-lived state is more stable than the original one, and the system as a whole (which 
has lost one state) is more stable than originally. The wavefunction of this long-lived state 
is mixed in those of the original states. 

The stabilization of the system due to Eqs. (Ell to fll9p occurs in an analog manner. 
Due to level repulsion, the two states separate from one another in energy, such that their 
interaction with one another is, eventually, of the same type as that with all the other distant 
states of the system. That means, each of the original states has “lost” its partner, also in 
this case; and the reproduction of the two originally neighbored states is prevented. As a 
result, the interaction of the states of the system via the environment is reduced (since all 
states are distant in energy), with the consequence that the system can be described well as 
a closed (almost stable) system. In difference to the case with width bifurcation, however, 
the number of states of the system as a whole remains unchanged. 

Eqs. (na to (I19p with purely imaginary and real coupling strength ca, respectively, will 
seldom be realized. They allow us however to receive analytical results (see Sect. ED and 
to understand the basic mechanism. Our numerical results for the more realistic cases with 
complex bj show the same effects. The interesting critical parameter range is that between the 
position of the EP and that of the maximum width bifurcation or level repulsion, respectively, 
as shown in the hgures with complex uj. In this parameter range, the wavefunctions are 
strongly mixed and varies suddenly from the value 0 at the EP, according to (E]), to the 
value almost 1, characteristic for almost orthogonal states at maximum width bifurcation 
and level repulsion, respectively. When ^ 1, the wavefunctions remain strongly mixed in 
relation to the original ones in all cases. 

The results are very robust and show the same characteristic features in all cases studied 
by us. They hold true for systems with loss (corresponding to decaying systems) and also 
for those in which gain may occur (by absorbing particles from the environment). They hold 
true also when gain and loss are balanced. 

It should be underlined here once more that a quantum system is really open and its 
properties are strongly influenced by the environment of scattering wavefunctions only in 
the vicinity of EPs. Here, neighboring states may strongly interact via the environment and 
may cause some decoupling of the whole system from the environment, as shown above. As 
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a result of this decoupling, the system is stabilized; behaves “linearly”; and Fermi’s golden 
rule is applicable. 

The eigenvalues shown in Fig. Elleft panel are real in the parameter range between the 
two EPs. This might be interpreted as a signature of PT-symmetry. The Hamiltonian 
is non-Hermitian also in this parameter range. The wavefunctions are biorthogonal 
and the phase rigidity is different from 1 for all parameter values including those for which 
the eigenvalues are real. The corresponding r* are near to 1, and not equal to 1. All 
calculations in this parameter range for realistic systems can therefore be performed, to a 
good approximation, by using a Hermitian Hamiltonian. Nevertheless, such a calculation for 
an open quantum system remains an approximation, although it will provide good results. 

Moreover, the parameters used in Fig. [3] left panel are unrealistic for a physical system. 
The coupling parameter u is usually complex as discussed in, e.g., {2, 13|. The eigenvalues 


obtained in a corresponding calculation are no longer real in a certain hnite parameter range, 
see the example Fig. |3]right panel. Also in this case, FT symmetry breaking may appear and 
the behavior of the system at maximum width bifurcation is determined by the nonlinear 
source term involved in the Schrodinger equation for an open quantum system; and we have 
the jump-like transition —)■ 1 at maximum width bifurcation also in this less symmetric 
case. 

The phenomenon of almost orthogonal wavefunctions at maximum width bifurcation (or 
maximum level repulsion) is robust as Figs. [Utoll] for different two-state systems show. It 
is not an artifact of the two-state model ([I]) since it appears also in calculations with more 
than two states (see Sect. IVBIl . 


B. N > 2 states 

All our calculations with N > 2 states are performed in the parameter region in which a 
higher-order EP is expected. Signatures of the existence of such a higher-order EP are not 
found in any of the results. This is, of course, not astonishing since every EP is a point in 
the continuum and therefore of measure zero. Also the second-order EPs (crossing points 
of two eigenvalue trajectories) can be identihed only by their influence onto observables in 
their neighborhood. Our results show clearly that this holds true also for higher-order EPs. 

One of the characteristic features of an EP (i.e. of the crossing point of two eigenvalue 
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trajectories) is that the two eigenfunctions in its surrounding are mixed due to the coupling 


of the system to a common environment 


13j. This environmentally-induced interaction of 


the states is large at and near to an EP (where the phase rigidity of the wavefunctions 
is reduced, as discussed in Sect. m- A nearby state does therefore not interact with the 
original states which cross at the EP. It interacts rather with some states, the wavefunctions 
of which are mixed in those of the original states (according to flTT]) j. At and near to these 
crossing points, new EPs of second order appear, the ranges of influence of which overlap. 
In other words, a clustering of EPs occurs, see the discussion in section HVl 

It is this phenomenon of clustering of EPs which we see in Figs. O to [71 It occurs in 
the parameter range of a higher-order EP. The results are generic and provide us valuable 
information on the dynamics of open quantum systems. Most interesting is the phenomenon 
of the jump-like enhancement of the phase rigidity when the maximum width bifurcation is 
parametrically approached. This effect is of the same type as that observed numerically for 
two states and discussed in detail in Sects. IV Al and IVI A[ respectively. 


VII. CONCLUSIONS 

In this paper we have described open quantum systems by means of a Schrodinger equa¬ 
tion the Hamiltonian "H of which is completely non-Hermitian. It contains explicitly (in the 
non-diagonal matrix elements) the interaction of the states via the environment. The eigen¬ 
values of "H are complex and the eigenfunctions are biorthogonal. Most interesting property 
is that the eigenvalues of two states may coalesce in one point (the so-called EP), at which 
also the corresponding eigenfunctions are the same, up to a phase a.TKeEP.ares.„g„.a. 
points and play an important role for the dynamics of open quantum systems. 

We used also the equivalent description of the system by means of a Schrodinger equa¬ 
tion with the non-Hermitian Hamilton operator "Ho (with vanishing non-diagonal matrix 
elements) and source term. Here, the interaction of the states via the environment is con¬ 
tained in the source term, and not in the Hamiltonian. The source term is nonlinear near 
and at EPs. It drives the behavior of the open quantum system and determines the dynamics 
of open quantum systems. 

Our main concern of the present paper is the phase rigidity r* of the eigenfunction of "H. 

This value provides a quantitative measure for the biorthogonality of the wavefunction of the 
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state i, i.e. for the possibility to influence the properties of the system by the environment. 
It holds 1 > Tj > 0. At Tj 1, the wavefunctions are almost orthogonal, very much like 
the eigenfunctions of a Hermitian operator. For vanishing however, the eigenfunctions of 
the non-Hermitian operator are really biorthogonal, and the influence of the environment is 
extremely large. This influence may cause, among others, a mixing of the wavefunctions of 
the different states via the environment. In [^, 13|, the modihcation of the eigenfunction <Fj 
of the non-Hermitian Hamilton operator "H due to its coupling to other states of the system 
via the environment is studied in detail. The resulting mixing of the wavefunctions can be 
expressed by the relation fITT]) . At and near to an EP, the mixing is extremely large. 

We studied hrst the phase rigidity in a two-level system around an EP. In all our calcula¬ 
tions, the coupling strength u between system and environment is hxed. Only the energies 
Ci or widths y* of the states are parametrically varied. We have ~ 1 far from an EP and 
Tj —)• 0 in approaching an EP. This result is expected from analytical studies. 

We observe however also another result in the critical region around an EP which is 
completely unexpected. In approaching the maximum width bifurcation and level repulsion, 
respectively, the value of the phase rigidity varies rapidly from its value r* < 1 to r* 1. 
That means, that the two wavefunctions are almost orthogonal when the width bifurcation 
or level repulsion is maximum. This jump-like variation of the phase rigidity is observed at 
hxed coupling strength between system and environment. It is caused therefore exclusively 
by the nonlinear source term of the Schrodinger equation. The wavefunctions of the states 
remain mixed at this critical parameter value, although they are almost orthogonal according 
to Ti —)■ 1. 

This phenomenon occurs not only in the simple two-state model (see Sect. IVAp but 
also in the case with more than two states (see Sect. IV Bj) . In the hrst case, we have 
well separated EPs while there is some clustering of EPs in the second case. That means, 
the clustering of EPs does not destroy the ehect, see Figs. [5] to [71 Quite the contrary, 
the clustering of many EPs causes a dynamical phase transition from an open quantum 
system (with biorthogonal eigenfunctions of its states) to an almost closed system (with 
almost orthogonal eigenfunctions of its states). The underlying process is irreversible and 
causes a stabilization of the whole system, meaning that the open system can be described 
approximately as a closed system. The wavefunctions at both sides of the dynamical phase 
transition are non-analytically related to one another and diher fundamentally from one 
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another. This feature is characteristic of any phase transition. The wavefunctions of the 
states on one side of the phase transition might be obtained by using the two-body residual 
forces derived from forces between free particles. This will be impossible, however, on the 
other side of the transition where the wavefunctions are modihed due to the mixing of the 
different states of the system via the common environment. 

Our results provide the following generic feature of open quantum systems. When two 
states are near to one another in energy or in lifetime, they may strongly interact with one 
another via the continuum of scattering wavefunctions due to the existence of a singular 
point (EP) in their vicinity. Here, the Schrodinger equation contains non-linear terms; 
irreversible processes occur; and the whole system will be stabilized. As a result, the system 
behaves very much like a closed system that is localized in space: The eigenstates are almost 
orthogonal and the eigenvalues are almost real (and sometimes even completely real j^). 
Such a situation can be described well by a Hermitian operator where the lifetime of a state 
does not appear explicitly. By this, the meaning of lifetime for the characterization of the 
individual states of the system is lost. Characteristic of the states are solely their energies 
and wavefunctions, while the lifetimes can be obtained by using perturbation methods. 

Nevertheless, the results presented in our paper show that energy and time are related 
to one another in quantum mechanical systems. It is shown in [15] that time is bounded 
from below in non-Hermitian quantum physics. This follows from the fact that the decay 
widths (inverse proportional to the lifetimes of the states) cannot increase limitless. Thus, 
time is bounded from below in the same manner as energy, in contrast to the assumptions of 
Hermitian quantum physics. Pauli has used this argument of Hermitian quantum physics, 
see e.g. [23], in order to conclude that the uncertainty relation between time and energy 


can, on principle, not be derived, for details see e.g. 


ITj . The uncertainty relation between 


energy and time remained therefore a puzzling phenomenon in Hermitian quantum physics. 
As our results show, this phenomenon is not at all puzzling in non-Hermitian quantum 
physics. 

n 

Concluding, we recall the phenomenon of resonance trapping [2[ observed many years 
ago. Resonance trapping occurring in an open quantum system coupled strongly to the 
environment, prevents the overlapping of individual resonance states. Consequently, the 
^stem is practically always in the regime of weakly (or not) overlapping resonances, see e.g. 


26l [ . In analogy to this phenomenon, an open quantum system can be described quite well 


26 






by a Hermitian Hamiltonian on both sides of the dynamical phase transition. This statement 
agrees completely with experience. Interesting non-trivial features of open quantum systems 
appear only in the parameter range in which a clustering of EPs and therewith a dynamical 
phase transition occurs. One of many examples is the relation between reduced phase rigidity 


and enhanced transmission through a quantum dot 

The results presented in this paper are generic. We believe that they will initialize 
further studies for concrete systems under concrete conditions. By this, they will provide 
new interesting results for open quantum systems, especially in the parameter range of a 
dynamical phase transition. 
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